Intermittency and the passive nature of the magnitude of the magnetic field 
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It is shown that the statistical properties of the magnitude of the magnetic field in turbulent 
electrically conducting media resemble, in the inertial range, those of passive scalars in fully de- 
veloped three-dimensional fluid turbulence. This conclusion, suggested by the data from Advanced 
Composition Explorer, is supported by a brief analysis of the appropriate magnetohydrodynamic 
equations. 
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In recent years, the problem of turbulent advection 
and diffusion of scalar and vector fields, both passive 
and active, has received renewed attention (see, for in- 
stance, Refs. I']-'^ and the papers cited there). Cor- 
responding experimental results are still few in number, 
partly because one is still grappling with the precise cor- 
respondence between the hypotheses of the theory and 
their experimental realization. We discuss here one of the 
most interesting examples considered recently: the case 
of magnetohydrodynamics (MHD) in which the magnetic 
field fluctuation B is described by the equation 
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V X (v X B) -f-r;V^B. 



(1) 



Here, v is the turbulence velocity and the rj is magnetic 
diffusivity. Equation (1) can be regarded as a vector 
analogue of the advection-diffusion equation 
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{v ■w)e + Dv^e 



(2) 



for the evolution of a passive scalar 6 subject to molecular 
diffusivity D. Aside from the fact that B is a vector and 
9 a scalar, the equations are different also because v in 
Eq. (1) can be affected quite readily by the feedback of 
the magnetic field B. Our interest here is to explore the 
extent of similarities, despite these obvious differences, in 
the inertial range statistics of the magnetic and passive 
scalar fields. 

Solar wind is an excellent natural "laboratory" for the 
MHD problem. It is known that the statistical properties 
of velocity fluctuations in the solar wind are remarkably 
similar to those observed in fluid turbulence 5] . It is also 
known that the plasma power spectra of the magnetic 
field and velocity fluctuations often contain an "inertial" 
range with a slope of approximately —5/3 (for reviews, 
see |6j,0). The approximately —5/3 power-law is es- 
pecially common for magnitude fluctuations B = \J Bf 
(the summation over repeated indexes is assumed) of the 
magnetic field, as one can see in Fig. 1. For comput- 
ing this spectrum, we have used the data obtained from 
Advanced Composition Explorer (ACE) satellite magne- 
tometers for the year 1998. In this period, the sun was 
quiet and the data are statistically stable. The range 
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FIG. 1: Energy spectrum of the magnitude B — yj Bf of the 
magnetic field B in the solar wind plasma, as measured by 
the ACE magnetometers in the nanoTesla range for the year 
1998 (4 min average). 



of scales for which the "—5/3" power holds is taken to 
be the inertial range; the smaller scales are obliterated 
because of the instrument resolution and the truncation 
at the large-scale end is governed by the record length 
chosen for Fourier transforming. 

The nature of the spectrum for each individual com- 
ponent of the magnetic field is more variable from one 
component of B to another, and from one situation to an- 
other, perhaps because of large anisotropics in the mag- 
netic field B, but the result for the magnitude of B seems 
more robust. Scaling spectrum with the —5/3 slope (Fig. 
1) is quite typical of that observed for passive scalar fluc- 
tuations in fully developed three-dimensional fluid tur- 
bulence (the so-called Corrsin-Obukhov spectrum |^). 
Spurred by this similarity, we were motivated to explore 
further the properties of the magnitude B, and compare 
them with those of the passive scalar. More detailed sta- 
tistical information is provided by the structure functions 
scaling 
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FIG. 2: Structure functions of magnetic field magnitude in the 
solar wind plasma as measured by the ACE magnetometers 
in nanoTesla for the year 1998 (4 min averages). The straight 
lines (the best fits) are drawn to indicate the scaling law (3) 
in the inertial range. 



FIG. 3: Scaling exponents (3) calculated for B in the solar 
wind (circles) and for the passive scalar in the atmospheric 
turbulence (triangles, Q), and in the direct numerical simu- 
lation of 3D fiuid turbulence (stars, 0) 



where 

ABr = B{t + T) - B{t). (4) 

The exponent (,2 is directly related to the spectral ex- 
ponent (in our case ^2 ~ 5/3 ~ I — 2/3 [3). If the 
dependence of Cp on p is nonlinear, it is well-known that 
one has to deal with intermittency. 

Figure 2 shows the scaling of structure functions for 
the solar wind data. Slopes of the straight-line fits in 
the apparently scaling region provide us the scaling ex- 
ponents Cpi these are shown in Fig. 3 as circles. Trian- 
gles in the figure indicate experimental values obtained 
for temperature fluctuations in the atmosphere 0. The 
other experimental data are in agreement with 

each other to better than 5%. The * symbols are for the 
passive scalar field obtained by numerically solving the 
advection-diffusion in three-dimensional turbulence 0. 
It is clear that the exponents for the passive scalar data 
are in essential agreement with those for the magnitude 
fluctuations of the magnetic field. 

One can, in fact, analyze the solar wind data somewhat 
differently using the notion of the extended self-similarity 
(ESS). Since, empirically, the fourth order exponent is 
quite closely equal to 1 for the magnetic field magnitude, 
i.e., 

{\ABr\')^r, (5) 

we can extend the scaling range (and consequently im- 
prove the confidence with which those exponents are de- 
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FIG. 4: Extended self-similarity (ESS) of the magnetic field 
magnitude in the solar wind plasma. The straight lines (the 
best fit) are drawn to indicate the ESS (6). 

termined) by redefining them as 

{\ABrn-{\ABr\')<^. (6) 

Figure 4 shows the ESS dependence (6). The slopes of 
the best-fit straight lines in this figure provide us with 
the ESS scaling exponents Cp, which are shown in Fig. 5 
as circles. The other symbols have remained unchanged 
from Fig. 3. The shift of the exponents in comparison 
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FIG. 5: The same as in Fig. 3 but using the ESS method (6) 
for B. 



If the statistical behaviors of 9 and B are to be similar, 
as suggested by Figs. 1, 3 and 5, we should be able to 
observe the underlying similarity between Eqs. (2) and 
(8). There is a major difference corresponding the pres- 
ence in Eq. (8) of the production term XB. However, 
given the empirical indications that B and 9 arc similar 
in the inertial range, it is appropriate to look for circum- 
stances under which the A term in Eq. (8) may be small. 
The second term in A is assured to be small because the 
smallness of the magnetic diffusivity 77, but difficulties 
may arise from the first term on the right hand side of 
Eq. (9). 

To eliminate directional dependencies in Eq. (8), let us 
make the following conditional average of that equation. 
That is, fix the magnitude B in the vector field B = Bn 
while performing the average over all realizations of the 
direction vector field n permitted by the vector equation 
(7). Let us denote this ensemble average as (...)„• From 
the definition, this averaging procedure does not affect 
B itself, but affects the velocity field v and the "friction- 
stretching" coefficient A in Eq. (8). We thus obtain 



to those from ordinary self-similarity is about 4%, but 
the scaling interval for ESS is considerably larger. This 
increased scaling range is well-known in other contexts 

The results shown in Figs. 3 and 5 suggest that at 
least up to the level of the fourth-order the scaling ex- 
ponents for the passive scalars and for the magnitude of 
the magnetic field are essentially the same. This is both 
surprising and thought-provoking, and needs to be un- 
derstood further. To this end, let us return to Eq. (1) 
and specialize jl3| . for simplicity, to the case of incom- 
pressible moving medium (V • v = 0). Equation (1) can 
then be rewritten as 
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'dt 



-(v- V)B + (B- V)v + ryV2B. 



(7) 



Let us now consider the equation for the magnitude B of 
the magnetic fluctuations given by B = Bn, where n is 
the unit vector with its direction along B: rii = Bi/B. 
Multiplying both sides of Eq. (7) by the vector n and 
taking into account that nf — 1 we obtain 
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'dt 



■V)B + riV^B + XB, 



(8) 



in which the "friction-stretching" (or the production) co- 
efficient A in the last term has the form 
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A = Uin-j—— 

OXi 



V 



drij 



(9) 



with the indices i and j representing the space coordi- 
nates, and the summation over repeated indexes is as- 
sumed. The first term on the right hand side of Eq. (9) 
is crucial for any dynamo effect. 
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T]V^B+{X)nB. (10) 



It is worth emphasizing that the solutions of the original 
equation (7) satisfy Eqs. (8) and (10), but not all possi- 
ble formal solutions of the Eqs. (8) and (10) satisfy Eq. 
(7); similarly, not all formal solutions of Eq. (10) sat- 
isfy Eq. (8) while all solutions of Eq. (8) do satisfy Eq. 
(10). Restricting comments to the relationship between 
Eqs. (8) and (10), the solutions of the two equations are 
the same only if the initial conditions are the same and 
if realizations of (v)n and of (A)n, related to these ini- 
tial conditions by the conditional average procedure, are 
obtained from solutions applicable to Eq. (8). 

Returning now to Eq. (10), the conditionally averaged 
velocity field (v)„ may posses statistical properties that 
are different from those of the original velocity field v, 
and there can be circumstances under which (A)„ — 0, 
or small. If so, the similarity between Eqs. (2) and (10) 
(and, consequently, Eq. (8)) can be the basis for the sim- 
ilarity in statistical properties of their solutions. There- 
fore, finding conditions under which (A)n = 0, or small, 
seems to be a useful exercise. 

It is, however, difficult to guess a prioiri when (A)n is 
negligible, because there is no small parameter for the 
stretching part of A. Therefore, let us consider a generic 
set of conditions, presumably for the inertial range, which 
can result in {nirijdvi/ dxj)ii = 0. This can be a combi- 
nation of isotropy, which yields 



and 



{n,nj)n = (z ^ j) 



{nl)n = {nl)n = {nl) 



(11) 
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and statistical independence 

{ninjip)n = {ninj)n{(p)n, (12) 

where ip = dv^/dxi for arbitrary k and 

We should emphasize that the conditional average in- 
dicated by (. . .)n and the global average indicated by 
(. . .) are quite different; because of this, the quantity B 
in (10) remains a fluctuating variable. To eliminate the 
stretching part from the conditionally averaged coeffi- 
cient (A)n — this being critical for explaining the observed 
similarity in the scaling of structure functions between B 
and 9 — one does not need to satisfy conditions (11) and 
(12) for all realizations of the magnetic field B, but only 
for the subset of realizations that gives the main sta- 
tistical contribution to the structure functions (3). Let 
us name this subset of realizations as /. The structure 
functions (3) depend on the statistical properties of the 
increments with respect to r, namely AS,-, belonging to 
the inertial range of scales. One of the consequences of 
intermittency is that the statistical properties of the in- 
crements are essentially different from those of the field B 
itself. Therefore, the subset /need not generally coincide 
with the subset G, say, that gives the main statistical 
contribution to the global average {niTijdvk/ dxi) . This 
means, in particular, that the conditions (11) and (12) 
can be valid for the inertial interval (i.e. for subset 7), 
while globally (i.e. for subset G) these conditions could 
well be violated. 

We now use conditions (11) and (12) in the presence of 
the incompressibility condition dvi/dxi = and obtain 

That is, the difference between the passive scalar equa- 
tion (2) and the conditionally averaged equation (10) for 
B is reduced to pure "friction" with the friction coeffi- 
cient given by (13). Equation (10) can then be reduced 
in Lagrangian variables to 

= (^)"^' (14) 

with the "multiplicative noise" (A)n given by Eq. (13). 
Weak diffusion of Lagrangian "particles" can be de- 
scribed as their wandering around the deterministic tra- 
jectories. Introduction of a weak diffusion is equivalent to 
introduction of additional averaging in Eq. (14) over ran- 
dom trajectories 14]. The small parameter ry in (13) and 
(14) will then determine a slow time in comparison with 
the time scales in the inertial interval and will therefore 



not affect scaling properties of B in the inertial interval. 
This explains the similarity of scaling between B and 6. 

In summary, we have shown that remarkable scaling 
similarities exist in the inertial range between the passive 
scalar and the magnitude of the magnetic field in MHD 
flows. Motivated by this observation, we have derived the 
dynamical equation for B and argued that, under circum- 
stances governed by Eqs. (11) and (12), dynamical simi- 
larity exists between the equations governing the passive 
scalar and the magnitude of the magnetic fluctuations. 
The conditions under which (11) and (12) are valid, per- 
haps best described as "directional randomness" , may 
be quite general and applicable to other circumstances 
(though, perhaps, not for the vorticity equation whose 
form is strongly nonlinear in any turbulent situation). 
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